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SYNOPSIS 


The  design  of  electron  guns  or  diodes  for  high  power  microwave  tubes  is  a  complicated  two 
dimensional  spatial  and  three  dimensional  velocity  space  problem.  It  involves  the  calculation  of  non¬ 
linear  particle  orbits  and  self  consistent  electric  and  magnetic  fields.  The  design  problem  is  expecially 
difficult  in  the  case  of  diodes  for  the  formation  of  helical  electron  beams  used  in  gyrotrons  because 
such  beams  must  have  high  transverse  momentum  as  well  as  low  velocity  spread  and  are  sensitive  to 
space  charge  forces.  Most  gyrotron  electron  gun  designs  are  based  on  the  use  of  electron  trajectory 
codes;  however,  the  complexity  of  the  problem  is  such  that  a  design  procedure  using  only  a  simulation 
code  can  be  very  tedious  and  expensive.  For  these  reasons  one  dimensional  approximations  to  the  flow 
are  extremely  useful.  Such  one  dimensional  flows  can  be  used  to  synthesize  first  approximations  to  the 
electrode  shapes.  The  purpose  of  this  paper  is  to  show  how  previous  work  on  this  technique  for  non- 
relativistic  beams  can  be  extended  to  the  relativistic  intense  electron  beam  regime. 

,  I 

The  analysis  is  based  on  four  one  dimensional  planar  or  conical  models  for  electron  flow  in 
diodes.  The  planar  models  are  relativistic  and  can  be  used  to  obtain  electrode  shapes  for  space  charge 
limited  flow.  The  conical  models  shed  light  on  important  geometric  effects.  Principally,  it  is  shown 
that  the  space  charge  limited  current  is  larger  in  realistic  geometries  than  in  the  corresponding  planar 
configuration.  A  design  procedure  based  on  these  flow  models  is  illustrated  by  application  to  a  diode 
for  a  1  GW  b^m  for  a  35  GHz  gyrotron. 
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ONE  DIMENSIONAL  MODELS  FOR  RELATIVISTIC 
ELECTRON  BEAM  DIODE  DESIGN 


I.  INTRODUCTION 

The  design  of  a  diode,  or  electron  gun,  for  a  microwave  tube  is  a  complicated  two  dimensional 
spatial  and  three  dimensional  velocity  space  problem.  It  involves  a  calculation  of  nonlinear  particle 
orbits  and  self  consistent  electric  and  magnetic  fields.  The  design  problem  is  especially  difficult  in  the 
case  of  diodes  for  the  formation  of  helical  electron  beams  used  in  gyrotrons  because  such  beams  must 
have  high  transverse  momentum  as  well  as  low  velocity  spread  and  are  sensitive  to  space  charge  forces. 
In  virtually  all  diode  design  work,  the  basic  tool  is  an  electron  trajectory  code  (Herrmannsfeldt  1979) 
which  solves  for  the  electron  orbits  and  the  self  consistent  electric  and  magnetic  fields  in  an  arbitrary 
two  dimensional  configuration.  Since  this  problem  is  so  complex  and  the  parameter  space  so  vast,  a 
design  procedure  using  only  a  flow  simulation  code  can  be  very  tedious  and  expensive.  Also,  since  the 
electron  trajectory  code  begins  with  vacuum  fields  as  a  first  approximation,  it  is  often  difficult  to  con¬ 
verge  on  flow  solutions  with  strong  space  charge  effects.  For  these  reasons,  one  dimensional  approxi¬ 
mations  to  the  physics  are  extremely  useful.  This  is  particularly  true  for  space  charge  limited  diodes, 
suitable  for  intense  relativistic  electron  beams  which  are  of  interest  tor  high  power  gyrotrons. 

The  basis  of  this  paper,  and  many  other  works  in  this  area,  is  that  there  are  useful  one  dimen¬ 
sional  approximations  to  the  orbit  in  an  infinite  medium.  Then  the  infinite  medium  approximation  is 
relaxed  by  taking  only  a  spatially  limited  part  of  the  infinite  flow  pattern,  and  focusing  electrodes  are 
used  to  create  the  fields  set  up  by  the  remainder  of  the  beam  which  was  excluded  (Marker  1960a, 
1960b;  Dryden  1962;  Kirstein,  Kino  and  Waters  1967;  Tsimiring  1972,  Manuilov  and  Tsimring  1979; 
and  Fliflet  et  al.  1982).  The  shapes  of  the  focusing  electrodes  are  determined  by  the  particle  orbit. 
Along  the  orbit,  both  the  potential  and  its  normal  derivative  are  known  Laplace's  equations  must  then 
be  solved  exterior  to  this  orbit  to  determine  the  equipotentials  (ihai  is.  the  shape  of  the  focusing  elec¬ 
trodes).  However,  the  solution  in  the  open  exterior  region  is  unstable  To  determine  the  solution,  a 
scheme  based  on  conformal  mapping  is  used.  The  orbit  equation  is  rewritten  as  an  equation  valid  in 
the  complex  plane.  Then  the  potential  is  calculated  by  for  real  r  and  imaginary  r  at  some  fixed  value 
of  real  r.  This  involves  the  solution  of  a  wave  equation,  which  has  stable  solutions  for  Cauchy 
boundary  conditions  on  an  open  surface.  Evaluating  the  solution  at  the  imaginary  part  ol  /  =  0  gives  the 
potential  as  a  function  of  r  at  the  fixed  value  of  the  real  part  ol  r  chosen.  Thus,  a  stable  solution  is 
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buill  up  by  taking  a  two  dimensional  projection  of  a  solution  in  three  dimensions.  The  equipotentials 
are  then  appropriate  places  for  the  location  of  focusing  electrodes. 

The  only  trouble  with  this  method  is  that  once  the  beam  becomes  relativistic,  and/or  self  mag¬ 
netic  fields  become  important,  as  is  the  case  with  intense  pulsed  beams,  there  is  no  one  dimensional 
model  which  describes  the  flow  This  paper  derives  a  variety  of  approximate  one  dimensional  models 
for  space  charge  limited  flow  in  a  diode  with  an  applied  magnetic  field.  While  no  model  is  exact  (the 
electron  flow  is  inherently  two  dimensional),  they  provide  important  insight  and  allow  for  approximate 
synthesis  of  the  electrodes.  The  actual  electrodes  can  then  be  quickly  perfected  with  the 
Herrmannsfeldt  code. 


Finally,  we  note  other  attempts  at  one  dimensional  models.  A  model  similar  to  one  of  those  in 
the  next  section,  but  with  the  magnetic  field  parallel  to  the  cathode  plane  has  been  derived  by  Gold¬ 
stein  (1976).  Another  approximau  model,  valid  in  planar,  cylindrical  or  spherical  geometry  is  the  Bril- 
louin  model  of  Creedon  (1975)  Here  the  self  magnetic  fields  are  assumed  so  strong  that  electrons 
flow  on  equipotential  surfaces,  so  that  the  model  does  not  describe  the  emission  from  an  electrode. 
Another  model  (Ott  et  al.  1977)  describes  the  electrostatic  field  in  two  dimensions,  but  the  electron 
flow  in  one  dimension.  Finally,  we  note  an  approximate  two  dimensional  analytic  model  of  focused 
electron  flow  in  diodes  (Goldstein  et  al.  1974). 

Section  II  describes  four  one  dimensional  models  for  electron  flow  in  diodes.  The  first  two  are 
planar,  relativistic  but  either  with  or  without  self  magnetic  fields.  The  second  two  are  conical,  with  no 
self  magnetic  fields,  and  either  nonrelativistic  or  superrelativistic.  The  nonrelativistic,  but  not  the 
superrelativistic  model,  is  well  known  and  has  been  used  in  several  applications.  The  last  model  can  in 
principle  accommodate  self  magnetic  fields,  although  we  have  not  exploited  this  feature.  We  rely 
mostly  on  the  planar  models  since  they  are  valid  over  the  full  energy  range.  However  there  are  still 
important  geometric  effects,  and  the  conical  models  shed  light  on  them.  Principally  it  is  shown  in  Sec¬ 
tion  III  that  the  space  charge  limited  current  is  larger  in  realistic  geometry  than  in  the  analogous  planar 
configuration.  Section  IV  reviews  the  electrode  synthesis  technique  and  derives  electrodes  for  a  rela¬ 
tivistic  beam  in  a  planar  configuration.  Finally,  Section  V  derives  a  final  electrode  configuration  using 
the  Herrmannsfeldt  code.  The  goal  is  to  achieve  a  high  quality  beam  for  use  in  a  gyrotron  oscillator  at 
35  GHz  and  with  power  in  the  hundred  megawatt  range.  Accordingly,  the  diode  has  a  magnetic  field  of 
1.8  -  2.0  kG  and  produces  a  600  kV,  2  kA  beam  with  o  =  v^/v,,  —  0.2.  In  the  diode  region  and  subse¬ 
quent  drift  space,  the  beam  has  virtually  no  variation  in  a.  To  achieve  the  higher  values  of  «,  the  mag¬ 
netic  field  must  be  compressed  to  about  25  kG.  In  doing  so,  some  variation  in  a  across  the  beam  is 
induced,  although  the  beam  quality  is  more  than  sufficient  for  efficient  oscillator  operation. 
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II.  ONE  DIMENSIONAL  MODELS 

In  this  section  we  describe  three  separate  one  dimensional  models  which  should  be  of  use  in  the 
design  of  a  wide  class  of  diodes.  These  models  are,  respectively,  a  relativistic  slab  model,  a  conical 
nonrelativistic  model,  and  a  conical  superrelativistic  model.  No  fully  relativistic  one  dimensional  coni¬ 
cal  model  was  found.  We  have  used  the  slab  relativistic  model  to  explore  a  large  area  of  parameter 
space  for  interesting  design  candidates.  This  model  is  used  in  the  electrode  synthesis  described  in  Sec¬ 
tion  IV.  The  conical  models  were  used  to  estimate  geometric  corrections  to  the  slab  designs.  st>ecifi- 
cally  the  current  required  to  ('btain  space  charge  limited  flow, 

A.  Slab  Relativistic  Model 

The  slab  diode  is  completely  described  by  three  constants  ol  motion  lot  c\ery  (larlicle.  namclc 
energy  and  the  canonical  momenta  and  /’.  Mere,  the  cathode  is  the  plane  \  *=  0  and  the  anode  i>  at 
A  =  L.  The  coordinate  z  corresponds  to  length  along  the  cathode  face  and  the  magnetic  field  H  i  .  in 
the  A  — z  plane.  [The  coordinate  c  corresponds  to  the  azimuthal  angle  </>  in  the  more  LonvciitiMiial 
cylindrically  symmetric  diode.]  For  any  particle  which  leaves  the  cathode  a  =•  0,  these  invaiianis  ol 


motion  are 


A  £  =  ( y  1 )  me-  -f-  V  <I>  ( .\ ' 


P,  -  y  mv, 


P.  ■=  y  mv .  * 


z.,l  -*-^,4,1x1 
c 


1 1  ■■  1 ^  (  x 


where  ,v'n.  Zn  give  the  initial  position  ot'  the  txarticle  I  he  consi.mi  oi  in"’|,,n  /'  in  ihc  i , iiiniiii  .il 


momentum  associated  with  i  in  the  gaiiec  wiih  ■(  =  t  'xi  ^  W  x 
constant  of  motion  /’.  is  the  x.'nnniv.il  iii.-'i’c  iituni  .inn,.,  ,.iicv!  x.io 


I  I  I  V  I  I  he  ihiid 

Ml  .nil  i|  hci  g.nigi  ,  V'  ii  1 1 


/<,  =  .4,(.vt  .4.  —  .-t.lx)  H,  1  Hci  lUsc  ol  ihf  g.iiigc  ih.mg"  ihi.  I’oi-son  hr.iiVci  .'iiw"cn 

F^  and  P.  is  not  zero  However,  ihi'  ■oimii  mi  !..•  n  .u'l'ln.i'i.  In  innh  g.uiccN  I  ind  i 

correspond  to  sell  lields  and  H,  .  H  .nt  Mu  cxitiiialM  .ippiml  l'^■ld''  In  I  qs  i  m  .m.l  ■ 

.q,  =  .4,  +  A,  .4.  =  .-i..  le  .  .X  rcprc-'  i  '-  i!i  i  ■  iK  'nii  /<  i  oiiNl.m'  .  h.r  e  tucii  I  '  I  n 


the  invariants  of  the  motion  aie  zero  on  tl  .  I'O  !•,■ 


o  ■!  'I'  (  I  .nc  iKo  /CIO  Iheic  Hci  ;uisi 


v.z  symmetry  only  one  orhit  needs  to  he  lo,-;.,;  '  r-  .:  anhom  Io-n  <  t  gciu  i  ilii\  it  i.m  h,i\c  i. 
From  Fqs.  (l)-(3)  and  the  definition  ol  >  vxc  on.l 


V  ,  I  V.  V  • 


>  I  X  M 


y(x)  mcy(x) 


Vj  U,^)  - 


y 

y(x) 


_ a _ 

mcy(x) 


(x). 


(6) 


where  11,  -  qBJmc  and  yCx)  -  1  -  q^/vax?^  from  (1).  Since  v,  is  chosen  to  be  positive  in  (4), 
these  equations  will  apply  only  up  to  the  first  turning  point  in  x.  This  is  exactly  the  condition  that  no 
orbits  cross. 


If  a  steady  source  of  current  jq  is  supplied  at  the  cathode,  an  equilibrium  state  will  be  set  up  with 
Jx  Jo  (by  V  •  j  -  0)  with  density 

_ 

■ 

Thus,  Poisson’s  and  Ampere’s  equations  give,  respectively. 


d^  4jryo 

dx^’  V,  ’ 

(7) 

d^Ai  iiTJo  V, 

(8) 

H 

> 

N 

d^Ay  4iryo  v^ 

(9) 

dx^  "  v,c 

dy/dx  -  v/v„ 

(10) 

dz/dx  —  v,/v,. 

(11) 

The  equations  (7)-(ll)  can  be  integrated  with  the  aid  of  Eqs.  (4)-(6)  using  x  as  the  independent  vari* 
able  as  an  initial  value  problem  from  x  -  0.  The  initial  conditions  imposed  are  y  «  0,  r  —  0, 

♦  -  0,Ay  »  /4,  -  0,  (all  by  convention),  -  0  (space  charge  limited  flow),  dAJdx  —  0  (fl^  —  0, 

corresponding  to  -  0  in  the  actual  diode)  and  dAJdx  -  Sjo  The  current  yo  can  be  adjusted  to  pro¬ 
vide  any  positive  potential  at  a  given  anode-cathode  gap  spacing,  ^(L)  -  V.  For  —  0  and  in  the 

nonrelativistic  limit,  these  equations  give  y  ~  z  —  A^  -  A^  —  0,  and  the  Langmuir-Child  solution 
***  ~  A  typical  numerical  solution  of  the  equation  for  nonzero  B^  is  shown  in  Fig.  1.  Time 

does  not  occur  explicitly  in  the  equations  we  integrate.  Nevertheless,  the  best  results  are  obtained  by 

using  a  uniform  time  step,  i.e.,  with  dx  =  v,rfr  This  gives  the  best  resolution  of  the  gyro-motion  in 

the  magnetic  field  and  the  singular  behavior  near  the  cathode,  where  the  solution  approximates 
Langmuir-Child,  even  with  B^  ^0. 


An  equilibrium  fully  symmetric  with  respect  to  y  and  z  is  obtained  if  /i,  is  a  constant  1  he  n.mi- 
thesis  discussed  in  Section  IV  provides  electrode  shapes  external  to  the  electron  beam  such  that  ihe 
beam  will  remain  symmetric  with  respect  to  y  and  z  if  Jo  is  uniform  in  a  finite  emitiiiii!  rcyion 
0  <  2  <  Lf  (but  -  oo  <  y  <  oo,  corresponding  to  0  <  ob  <  27r).  However,  it  is  not  i)ussible  to 
obtain  the  correct  A,  by  such  a  method,  because  B,  is  proportional  to  the  curieni  /.  withm  ihc  oil'ii 
(see  Fig.  1.).  For  v_,  >  0  this  is  zero  for  the  bottom  orbit  (from  the  right  end  ol  the  emitting  region » 
and  increases  for  orbits  emitted  further  to  the  left,  but  can  nc'.er  be  equal  to  ihe  \alue  obtaiiKd  m  ihs 
pure  slab  limit  for  z  >  L<,,  because  some  current  is  thrown  out.  .Mternatively.  we  could  note  that  w  hh 
a  finite  emitting  region.  By  becomes  a  function  of  .v  and  In  older  to  bracks;  'h-  '  -i  t  soluium  lor  i 

planar  diode  with  a  finite  emitting  region,  we  have  iniegiatcd  the  equations  i  '  i  >  with  the  tun  i  '  .  ' 
and  again  with  -4;(a)  =  0.  The  first  model  is  appropiiate  when  the  emittinv  of  le.igiii  1  ■ 

long  or  for  orbits  near  the  top  of  the  beam  before  they  have  traveled  a  tlis,.iii>.e  Ar  equal  to  1  ,  i 
second  model  is  appropriate  for  L,.  — *  0  or  for  orbits  near  the  boiiom  ol  the  beam  RcmiIis 
from  these  models  will  be  discussed  in  Sec.  Ill.  The  effect  of  a  finite  emitting  region  on  t  .  A  .md 
for  that  matter,  the  effect  of  J,  upon  B.  for  the  L,.  —  oo  case,  is  negligible  lor  the  designs  we  siud'> 
because  of  a  strong  applied  B.. 


B.  Conical  Nonrelativistic  Model 

A  scaling  that  reduces  the  nonrelativistic  conical  diode  problem  to  one  dimension  is  well  kno.in 
but  we  include  a  brief  presentaiion  for  completeness  We  use  spheric. i'  .  s.idman  s  (rw  .,6/  i  .t;  .  ..ne 
a  diode  with  cathode  at  0-0,..  anode  at  «  -  and  svmme'i'.  fhe  lelc'.ini  cqu.iUi  O'  o  , 

B;  —  Bo  =  const,  and  =  0  (i  e  ,  ignorm;  self  an  |  /  i  are  . . . .  > ' 

”1  *  ■  • 'I  I  " '  o  '  I . ' 
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(16: 
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■  ^1  IJI  P  I  ^  IJI  I  I  I  ^1  I  BM  I  I  ■! 
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“T"^  (/"^wv,)  +  - : -  — 

r^  dr  r  sin  9  d9 


(sin  9  n  \g)  •’  0. 


These  equations  can  be  reduced  to  a  one-dimensional  form,  with  9  the  independent  variable,  by  assum¬ 
ing  the  scaling 


«!>  -  4>o  (f ),  (17a) 

V  -  r  Vo  (9),  (17b) 

n  -  no  (d).  (17(.) 

Equations  (13),  (14)  then  reduce  to 

y  m  Vo^  (9)  -I-  1?  <I>0  (6>)  =  0,  (18) 

v»o^/(?/^/«-(-2v,ol2(0)  =  O,  (19) 


where 


Q(9)  =  sin  9  v^o  +  sin^  0/2 
and  fto  ”  BqJ me.  Equation  (14)  gives 


‘’M  (»)  +  VrO  -  Vg^o  -  vio 


± 

m 


200  (0)  + 


v^o  (®)  Bq  sin  9 


From  (15)  and  (16)  we  obtain 


1 


sin  9  do 


sin  9  +  6  Oo  (fl)  =  —  ^irq  rig  (9), 


I 


sin  9  do 


sin  9  no  (0)  V(,o(@)j  -I-  3  no(9)  v,o  (9)  =  0. 


(20) 


(21) 


(22) 


(23) 


We  integrate  (19),  (21),  and  (22)  with  respect  to  9  to  obtain  v^o(9),  and  Oo,  respectively,  in  addi¬ 
tion  to  finding  r(9),  </>(9)  by 


do 


(Inr) 


v.o 

V(.0  ‘ 


dtj)  _  '^ito 

do  sin  v„o 


(24) 

(25) 


Finally,  no  and  v,o  are  obtained  from  (23)  and  (18) 


Note  that  self  magnetic  fields  cannot  be  incorporated  min  these  sin'iliirii\  soIuiiotis  'in.  r  o;. 
dynamical  eciuations  (13),  (141  require  B  =  B,,  (#<),  j  nq  \  =  r  j.  to)  fi.c  .  ( |  7h»  and  ( I  7u  >  I  v  i  . 
cannot  satisfy  V  x  B  =  4;?  j/c 


V. 


C.  Conical  Superielatisistic  Limit 


In  the  limit  ut  supeiielaiiMstic  electron  eneigies.  v  >  •  1  I  us  i  i  "’i-i  l-i!  .nc  u  pl.ii.,  d  I  \ 

y  mc'^  +  q  0 


d  •  d 

V ,  —  4  •  - 

r  sm  0  /!,  •  -r  q  , 

dr  r  dd 
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where  p=  >■  m  \  end  %  =  (I  +  p',  m  ’  hemmev  ■>  ---  'p: me.  i  ,  Im  l  ijiieii'  .  'i 

are  unchtinged  ( )ne  dimerision.d  egudiions  m  o  are  oht.iiiicd  in  this  liniii  h>  the  (ollowntg  ■  de 


=  r  4>|,  to) 
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P  =  r  Po  (d) 
V  “  Vfl  (0) 


V  r  >„(«• 

B  =  B.)  (»)  ___  '  d 

A  -■  r  A(,  lot  '■.! 

n  =  tin  t0)/r.  ■  .  u, 


Note  that  ri’-  end  ''7  x  B  both  scale  as  1/r.  so  that  self  fields  can  be  included  in  these  sohitirins  I  .|i!  ■ 
tions  (2f>),  1 2')  give 

yo'd*  nie’  e  q  <J>„  to)  =  0,  : 

Pun  ^PrnQ  0-  ■ 


Equations  (28),  (15)  and  (16)  give 


dPrO  , 

''tfO  '^rO  PrO  ^«QPeO  P*0 


(33) 


^  -n  h  *  ^ 

^^(sii,«»,V„)  +  „,v,„-0.  (35) 

The  9  component  of  Ampere’s  law  gives  B^o  —  —  no9  and  the  r  component 

^  (“"  *  «.»|  -  «  vk, 

is  consistent  by  (35).  The  <f>  component  of  Ampere’s  law  gives 

Bgo  “  ~  2  (36a) 

'de  ^  -  »o  9  v^o-  (36b) 

Equations  (30)-(36),  together  with  (24)  and  (25)  can  be  integrated  as  initial  value  equations  from 
0  »  00  as  in  the  planar  and  nonrelativistic  conical  diode  cases.  When  the  self  magnetic  fields  are  negli¬ 
gible,  B^,q  0,  5*0  “  ~  2  ^40  =  -  Bq  sin  0  and  (36b)  is  not  used. 

III.  APPLICATIONS  OF  ONE  DIMENSIONAL  MODELS 


In  this  section  we  describe  how  the  one  dimensional  models  of  Sec.  II  are  used  in  diode  design. 
The  slab  relativistic  model  has  been  used  primarily  to  scan  the  parameter  space  for  reasonable  design 
candidates.  Electrode  synthesis  has  been  applied  to  the  promising  designs,  with  further  electrode  shap¬ 
ing  required  both  for  final  tuning  of  the  design  and  for  other  practical  considerations.  This  model  has 
been  used  to  shed  light  on  certain  aspects  of  the  design,  for  example,  the  dependence  of  the  results  on 
the  self  B^.  The  conical  models  have  been  used  to  clarify  the  geometric  effects,  specifically  the  varia¬ 
tion  of  the  space  charge  limited  current  obtained  by  the  Hermansfeldt  code  from  the  value  predicted  by 
the  slab  model.  Let  us  assume  that  we  want  to  produce  a  beam  with  a  =  \pil P\\  \  =1  in  an  interaction 
region  of  radius  ry  -  1.4cm,  with  a  field  B.j  -  24  kg.  Assuming  that  the  beam  is  adiabatic  in  the  drift 
region  between  the  diode  and  the  interaction  region,  i.e.. 


B, 


(37) 


and  taking  a,  —  0.2,  we  obtain  fl/  =  1.85kg,  r,  —  5.05  cm. 

We  have  used  the  slab  relativistic  model  to  compute  a  at  the  anode  as  a  function  of  the  angle  xo 
between  the  normal  to  the  cathode  and  the  externally  imposed  magnetic  field,  i.e., 

The  gap  voltage  was  <I>  -  600  kV,  and  the  current  density  Jo  was  70  A/cm^.  The  results, 
with  and  without  self  Bj,(x),  are  shown  in  Fig.  3.  Without  B^,  xo  ”  30°  gives  a,  =  0.2,  whereas  with 
By,  40°  is  correct.  This  indicates  that  self  By  plays  an  important  role  in  determining  a,  although  it 
has  little  influence  on  the  gap  spacing  d.  Since  By  (or  B^)  is  zero  for  the  bottom  orbit,  and  since  By 
(B^)  can  be  near  its  limiting  one  dimensional  value  for  orbits  near  the  top,  this  indicates  that  electrode 
synthesis  may  not  by  itself  produce  a  uniform  a  across  the  beam.  We  will  return  to  this  issue  in  Sec. 
V.  We  also  show,  in  Fig.  4,  the  dependence  of  a  and  d  upon  Jq,  for  xo  “  40°,  d)  =  600  kV,  including 
self  By(x).  It  is  clear  that  both  of  these  quantities  depend  critically  upon  Jq. 

As  we  shall  discuss  further,  the  Hermannsfeldt  code  in  this  parameter  range  shows  space  charge 
limited  flow  at  20-50%  higher  current  than  that  indicated  by  the  slab  model.  In  order  to  understand 
this  discrepancy  and  have  more  confidence  in  our  results,  we  have  investigated  this  effect  with  the  non- 
relativistic  and  superrelativistic  conical  models.  We  use  these  models  in  the  following  manner:  we  fix 
the  cathode  angle  flo  (see  Fig.  2),  a  potential  a  gap  spacing  d,  a  field  B^^,  and  a  radius  R  where  a  ray 
is  to  be  emitted.  We  adjust  Jq  until  the  potential  equals  d>  when  the  gap  spacing  'x%  d.  Then  we  vary  R. 
The  results  for  different  R  values  are  not  obtainable  from  the  scalings  (17)  or  (29)  because  we  do  not 
allow  <I>  and  d  to  scale  appropriately  with  /?  (<I>  ~  R^,d  ~  R  from  (17)or<I>~/?,</~/?  from  (29)j. 
For  a  case  with  d)  -  1.07  MV,  B^a  -=  2.4kg,  d  —  3cm,  “  )58  °,  we  obtain  the  results  shown  in  Fig.  5. 
For  both  the  nonrelativistic  and  superrelativistic  models,  the  results  fit  curves  of  the  form  Jq  =  y'o 
-^AfR.  and  the  limiting  value  agrees  with  corresponding  slab  nonrelativistic  or  slab  superrela¬ 

tivistic  models.  For  the  cases  considered,  Jq  for  /?  —  5  cm  is  20  to  50%  higher  than  for  the  slab  model, 
in  agreement  with  the  results  obtained  using  the  Hermannsfeldt  code. 

IV.  SYNTHESIS  TECHNIQUE  FOR  DETERMINING  ELECTRODE  SHAPES 

The  electrode  synthesis  technique  is  a  method  of  calculating  electrode  shapes  which  provide  lam¬ 
inar  flow  for  a  beam  with  self  electric  fields.  The  method  was  initially  formulated  by  Marker  for  planar 
and  axially  symmetric  cases  of  space  charge  limited  nonrelativistic  flow.  The  method  was  extended  to 
temperature-limited  MIG  type  guns  by  Manuilov  and  Tsimring.  The  general  approach  involves  finding 
a  set  of  ordinary  differential  equations  to  represent  the  beam  flow.  These  equations  are  used  to  obtain 
the  boundary  conditions  for  integrating  Laplace’s  equations  in  the  region  outside  the  beam.  A  difficulty 


in  the  direct  inplementation  of  this  procedure  is  that  the  present  problem  involves  Cauchy  boundary 
conditions  on  an  open  surface  and  the  solutions  of  Laplace's  equation  are  unstable  for  these  conditions. 
An  elegant  method  has  been  developed  by  Marker  which  reformulates  the  mathematical  problem  in  a 
way  which  yields  stable  numerical  solutions.  This  section  outlines  the  synthesis  technique  for  planar 
geometry  and  discusses  an  application  for  the  case  of  planar  relativistic  flow. 


The  synthesis  problem  consists  of  finding  the  electrostatic  potential  distribution  in  a  region  exter¬ 
nal  to  a  finite  laminar  flow  beam.  The  potential  satisfies  Laplace's  equation  outside  the  beam  or 
equivalently,  the  electrostatic  field  satisfies  the  Maxwell  divergence  and  curl  equations 


dE,  dE, 
— ^  +  — - 
dx  dz 


0 


(39) 


dE,  bE, 
bz  bx 


(40) 


with  the  boundary  conditions 

E,  -  =  E„  (41) 

on  the  outermost  trajectory  of  the  beam. 


Equations  (39)  and  (40)  form  a  system  of  elliptic  partial  differential  equations  whose  solutions  are 
unstable  for  the  present  case  of  an  open  surface  and  Cauchy  boundary  conditions.  However,  a  stable 
solution  can  be  obtained  by  the  following  approach  due  to  Marker. 


First,  the  beam  edge  trajectory  (x,  -  x,(/),  z,  -  z,(r)),  is  converted  into  the  straight  line  «  =  0 
of  the  (t.u)  plane  by  means  of  the  conformal  transformation 


z  +  ix  ~  Zf  if  +  iu)  +  iXf(f  +  iu). 


(42) 


This  can  be  done  because  the  equations  for  the  beam  trajectory  can  be  expressed  in  the  complex  plane 
by  analytic  continuation.  Since  coordinates  related  by  a  conformal  transformation  satisfy  the  Cauchy- 
Riemann  conditions. 


bx  _  bz 
bu  bt 


(43) 


bz  ^  ^ 

bu  bt ' 


Equations  (39)  and  (40)  can  be  expressed  in  the  form 


(44) 


•:v:;:> 

-«*Vi 


8£,  a£, 

du 


These  equations,  together  with 


d<t‘ 

du 


dt  ■ 


—  £  —  +  £ 
*6/  '  dt 


(46) 


(47) 


enable  calculation  of  the  potential  4>  in  the  transformed  plane  (r.u). 


The  second  step  is  the  transformation  of  the  elliptic  system,  Eqs.  (45)  and  (46)  in  the 
plane,  by  means  of  the  analytic  continuation 


‘S*v 

.vy.* 

*  N^V 

V"  ■ 


t  —  p  +  /q 

For  fixed  p  this  leads  to  the  hyperbolic  system  in  the  (q,«)  plane: 

d£^  _  .8£. 
du  ‘  dq 

d£j  i  dEx 

du  dq 

du  dq 

di  ^  idx 
du  dq 

Equations  (49)  and  (52)  have  a  stable  solution  for  the  present  boundary  conditions. 


(48) 

(49) 

(50) 

(51) 

(52) 


The  procedure  for  obtaining  the  equipotentials  is  illustrated  in  Fig.  6.  It  involves  solving  the  sys¬ 
tem  (49)  and  (52)  in  triangular  regions  such  as  ABC  in  Fig.  6.  To  obtain  the  solution  for  this  region  it 
is  sufficient  to  specify  Cauchy  boundary  conditions  on  the  line  AB.  These  conditions  are  the  analytic 
continuation  of  Eq.  (41)  and  are  formed  by  integrating  the  flow  equations  along  the  real  axis  from 
f  =•  0  to  /  -  £o  and  then  along  the  line  AB  by  means  of  the  substitution 


dt 


(53) 


'•'.'•■‘'I' 


Solving  Eqs.  (49)-(52)  in  the  region  ABC  by  a  finite  difference  method  allows  the  potential  to  be 
obtained  on  the  line  CD,  the  only  region  of  physical  significance.  By  translating  the  triangle  ABC  to 
other  values  of  p,  the  potential  distribution  over  the  entire  single  valued  region  of  the  plane 
(p,u)  -  (f,«)  can  be  found.  Applying  the  transformation  (42)  yields  the  equipotential  surfaces  as  a 
function  of  x  and  z. 


r-' 


II 


Electrodes  shapes  calculated  by  the  synthesis  technique  are  shown  in  Fig.  7  for  the  planar  rela¬ 
tivistic  flow  model  discussed  in  section  II  A.  The  synthesis  calculation  was  based  on  a  current  density 
of  70  A/cm^,  an  angle  of  40°  between  the  cathode  normal  and  the  external  magnetic  field,  and  an 
external  magnetic  field  by  1.85  kG.  The  anode-cathode  gap  voltage  is  600  kV.  Upper  beam  edge  tra¬ 
jectories  are  shown  with  and  without  the  approximation  self  magnetic  field  effect  included.  As  shown 
the  effect  is  quite  small  for  the  present  parameters  and  there  is  negligible  effect  on  the  calculated  elec¬ 
trodes. 

V.  FINAL  DIODE  DESIGN 

The  final  design  for  the  diode,  in  several  configurations,  was  achieved  by  means  of  an  electron 
trajectory  (Hermannsfeldt)  code.  The  electrodes  found  by  synthesis  generally  provide  a  beam  in  which 
a  varies  by  ±  50%  from  top  to  bottom.  This  variation,  which  is  not  present  in  the  slab  model,  is  due 
to  geometric  effects.  One  of  these  effects  is  the  fact  that  the  bottom  orbit  has  =  0,  whereas  the  top 
orbit  has  a  value  of  which  can  be  nearly  equal  to  the  one-dimensional  value  (for  fat  beams).  This 
effect  has  been  studied  in  Sec.  Ill  and  found  to  provide  just  such  a  variation  in  a. 

Our  first  design  is  for  a  600  kV,  2  kA  gun  which  can  use  a  minimal  amount  of  focusing,  depend¬ 
ing  upon  intercepting  the  outer  two  thirds  of  the  beam  at  the  anode.  This  extra  charge  (and  current) 
takes  the  place  of  some  of  the  focusing.  A  successful  design  of  this  type  is  shown  in  Fig.  8.  The 
cathode  face  is  at  40°,  which  we  found  in  Sec.  Ill  gives  a  =  0.2  with  =  1-85  kg.  The  current  den¬ 
sity  7o  =  70A/cm^  in  the  slab  model  gives,  for  a  cathode  surface  of  radius  5  cm  and  emitting  length  2.1 
cm,  4.6  kA.  Because  of  geometric  effects  as  discussed  in  Sec.  Ill,  we  used  6.0  kA,  corresponding  to  a 
perveance  k  =  13  micropervs.  This  is  near  the  space  charge  limit;  for  Ar  >15,  the  results  begin  to 
show  serious  signs  of  lack  of  convergence  usually  associated  with  approaching  the  space  charge  limit  at 
some  point  on  the  cathode  face.  The  results  show  a  very  flat  a  =  0.2,  ±10%  and  very  little  evidence 
of  orbit  crossing.  It  is  possible  to  reshape  the  electrodes  to  have  less  focussing  in  order  to  decrease  the 
electric  fields  on  the  parts  of  the  cathode  where  we  wish  to  inhibit  emission.  In  that  case  a  would  not 
be  as  flat  across  the  whole  6  kA  beam,  but  that  is  of  no  consequence.  However,  it  appears  that  the 
electric  fields  in  the  design  in  Fig.  6  are  below  300  kV/cm,  and  the  surfaces  of  anodized  aluminum  on 
the  focusing  electrodes  should  inhibit  emission  for  at  least  50  ns. 

In  Fig.  9  we  show  a  design  with  an  aperture  in  the  anode  allowing  roughly  the  middle  third  of  the 
beam,  and  a  short  drift  region  bounded  by  anode  surfaces.  There  is  very  little  orbit  crossing  still,  and  a 
is  quite  flat  in  the  central  third  of  the  beam.  (For  rays  intercepted  by  the  anode,  the  value  shown  for  o 
is  the  value  at  interception.)  However,  o  is  considerably  higher  now.  in  the  range  0.36  <  a  <  0.41. 
This  is  apparently  due  to  a  combination  of  effects,  including  the  space  charge  of  the  beam  and  the  finite 
Larmor  radius  of  the  orbits. 


An  alternate  design  for  a  600  kV,  2  kA  diode  is  shown  in  Fig.  10.  Here,  uniformity  of  the  beam 

is  achieved  by  focusing  electrodes  alone;  all  of  the  beam  is  allowed  to  pass  through  the  aperture  in  the 

anode.  The  electrode  shapes  are  nearly  identical  to  those  of  the  design  in  Figs.  8  and  9,  except  that  the 

emitting  region  is  one  third  the  length,  or  7mm;  the  perveance  is,  correspondingly,  k  =  4.3/x  micro- 

perus  and  the  magnetic  field  is  a  nearly  uniform  1.85  kg.  Again,  the  electrode  shapes  are  somewhat 

different  from  those  obtained  by  synthesis,  for  geometric  reasons,  and  the  maximum  electric  field  is  of 

order  300  kV/cm.  Figure  10  shows  essentially  no  orbit  crossing  in  the  diode  region  and  a  very  uniform 
a  “  0.19±5%  across  the  beam  at  the  anode  and  at  the  end  of  the  drift  region.  A  combination  of  coils 

produce  a  quite  uniform  2.2  kg  magnetic  field  in  the  diode  region  that  increases  to  20  kg  at  the  far  end 

of  the  compression  region.  The  average  value  of  a  at  the  far  end  of  the  compression  region  agrees  well 

with  the  adiabatic  value  —  1.2  based  on  a,  =  0.22,  B.  =  2.2  kg  at  the  beginning  of  the  compression 

region  and  5/  =  28kg  at  the  end.  However,  there  is  much  more  variation  in  a  across  the  beam  than 

adiabatic  theory  predicts;  from  (37)  one  can  show 

In  a  r  2 

-J, - -  =  . 

d\n  1  +  af 

which  shows  that  the  relative  variation  in  should  only  be  twice  the  relative  variation  in  a,  for 
■X,  <  <1.  This  is  a  common  effect  in  this  type  of  simulation  and  is  apparently  due  to  space  charge 
effects  that  become  more  pronounced  as  the  beam  slows  up.  The  results  shown  in  Fig.  10  a,c,d,e  are 
with  25  rays;  results  with  15  rays  as  shown  in  Fig.  10b  give  very  similar  results,  including  the  variation 
of  a  across  the  beam. 

Tests  have  been  made  to  determine  the  sensitivity  of  the  diode  performance  to  variations  in  B. 
and  the  potential.  For  these  tests,  a  600  kV,  2  kA  diode  with  B.  —  1.85  kg  was  used.  These  tests  are 
important  in  order  to  be  able  to  tune  an  actual  device.  In  addition,  it  is  important  to  know  whether 
variation  of  the  voltage  during  the  pulse  will  have  deleterious  effects.  The  nominal  diode  design  here 
has  a  =  0.25  across  the  beam.  Dropping  the  potential  to  450  kV  causes  the  beam  to  enter  the  aperture 
in  the  anode  less  than  a  centimeter  below  the  point  where  a  600  kV  beam  enters,  so  that  it  appears  that 
a  pulse  with  450  kV  <  <l)<  600  kV  will  produce  a  beam  that  can  still  fit  through  the  anode  aperture. 
Furthermore,  a  is  nearly  unchanged.  These  results  are  not  sensitive  to  the  perveance.  For  B.  ten  per 
cent  below  the  nominal  value,  i.e.,  1.67  kg,  the  beam  begins  to  scrape  the  top  of  the  drift  cavity,  and  a 
is  larger,  in  the  range  0.32  <  a  <  0.38.  For  B^  twenty  percent  above,  i.e.,  B.  =  2  22  kg  has  a  =  .20 
Thus  it  appears  that  a  scales  as  1/fl.,  so  that  varying  B.  may  be  an  effective  way  to  obtain  a  desired  a 
in  this  type  of  diode. 

An  important  consideration  for  realizing  the  design  based  on  focusing  electrodes  involves  the  con 
trol  of  emitting  regions  under  high  voltage  conditions.  If  the  electric  field  is  too  high,  everything  will 


ultimately  emit,  including  the  focusing  electrodes.  The  key  is  to  keep  the  field  sufficiently  low  and  the 
emissivity  of  the  emitter  sufficiently  high.  For  instance,  it  has  been  determined  (Kirkpatrick  et  al., 
1984)  that  with  focusing  electrodes  of  anodized  aluminium  and  emitting  surfaces  of  reactor  graphite, 
currents  in  the  kiloampere  range  could  be  generated  for  30  nsec  without  emission,  and  that  the  focus¬ 
ing  electrodes  do  not  emit  for  fields  as  high  as  600  kV/cm. 
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Fig.  1  —  Geometry  for  the  slab  relativiatic  diode  model.  For  the  bottom  orbit,  is  zero. 
For  the  top  orbit,  B^  is  the  value  obtained  by  the  slab  model  for  electrons  to  the  left  of  the 
dashed  line.  To  the  right,  B,  is  leal  than  the  slab  model  value. 
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Fig.  2  —  Geometry  of  conical  diode.  For  a  finite  emitting  region,  the  same  comments  hold 
for  self  B^.  The  angle  of  Fig.  I  corresponds  to  0g  -  n/2. 
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Fig.  3  —  Dependence  of  a  -  \pjpn\  and  (he  gap  spacing  d  upon  xo  slab  relativistic  model 
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Fig.  8  —  A  600  kV,  2  kA  diode  design  with  the  outer  two  thirds  of  s  6  kA  beam  to  be  intercepted  by 
the  anode.  In  (a)  are  equipotentials  and  actual  trajectories  integrated  by  the  code;  in  (b)  is  shown  a 
vs  ray  number,  from  the  bottom  oE  the  beam  to  the  top,  for  the  electrons  when  they  hit  the  anode. 
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Fig.  10  —  (■)  Diode  design  with  600  kV,  2  kA,  where  the  beam  uniformity  is  provided  by  focusing 
electrodes  alone,  i.e.  no  interception  of  beam;  (b)  a  at  the  aperture  in  the  anode;  (c)  a  at  the  end 
of  the  drift  region;  (d)  orbits  in  the  compression  region,  (e)  a  at  the  end  of  the  compression 
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Fi|.  10  (Continued)  —  (a)  Diode  design  with  600  kV,  2  kA,  where  the  beam  uniformity  is  provided  by  focusing  elec¬ 
trodes  alone,  i.e.  no  interception  of  beam;  (b)  a  at  the  aperture  in  the  anode;  (c)  a  at  the  end  of  the  drift  region;  <d) 
orbits  in  the  compression  region,  (e)  a  at  the  end  of  the  compression  region 


